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Abstract 

A search for single top quark production via flavour changing neutral currents (FCNC) was 
performed with data collected by the OPAL detector at the e^e~ collider LEP. Approxi- 
mately 600 pb"^ of data collected at y/s = 189-209 GeV were used to search for the FCNC 
process e+e~ —>■ tc(u) bWc(u). This analysis is sensitive to the leptonic and the hadronic 
decay modes of the W boson. No evidence for a FCNC process is observed. Upper limits 
at the 95% confidence level on the single top production cross-section as a function of the 
centre-of-mass energy are derived. Limits on the anomalous coupling parameters and kz 
are determined from these results. 
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1 Introduction 



In the mid 1990's, the LEP colhder at CERN entered a new phase of operation, LEP2, 
with the first e"'"e~ colhsions above the W"'"W~ threshold. Between 1998 and 2000, with the 
installation of additional super-conducting radio-frequency accelerating cavities, the centre- 
of-mass energy of the LEP collider was further increased. The LEP2 data accumulated 
at centre-of-mass energies between 189 GeV and 209 GeV have opened up a new kinematic 
domain for particle searches. 

The top quark mass was measured at the Tevatron collider to be 174.3 ± 5. 1 GeV/ [|I|, |^ . 
Due to this high mass, top quarks may only be singly produced at LEP2. Single top quark 
production in the Standard Model (SM) process e"''e~ — > e~t'etb has a cross-section of about 
10~^fb at LEP2 energies and can not be seen with the available luminosities. Another 
possible process for single top quark production is the flavour changing neutral current 
(FCNC) reactionQ: 

e+e-^tc(u). (1) 



Such FCNC are known to be absent at the tree level in the SM but can naturally appear 
at the one-loop level due to CKM mixing which leads to cross-sections of the order of 10~^ fb 
at LEP2 energies 0. Extensions of the SM such as supersymmetry, exotic quarks, and 
multi-Higgs doublet models could lead to an enhancement of such transitions H, 0, §]. 
In this paper the search for single top production via the FCNC reaction e"'"e~ — > tc(u) is 
reported. 

At the Tevatron, the CDF Collaboration performed a search for FCNC in the top de- 
cays t 7c(u) and t Z c(u) in pp collisions at a centre-of-mass energy of 1.8 TeV. They 
obtained upper limits at the 95% confidence level (CL) on the branching fractions 0: 
Br(t ^ C7) + Br(t U7) < 3.2% and Br(t ^ cZ) + Br(t -> uZ) < 33%. 

The FCNC reaction can be described with the parameters and kz which represent 
the tree-level 7 and Z exchange contributions to e^e~ — tc(u). Thus, the Born-level cross- 
section for single top production in e+e~ collisions for ^/s > rrit can be written as [^: 



cr e e 



tc(u)] 



rriT 



+ 1 h Sk^Kz 

4sin^ 2^9w(l-^)' 



2 2_f_ 

'^7^<? 2 



1 + 



2ml 



sin2 2^w(l-^)J 



(2) 



where s is the centre-of-mass energy squared, a is the fine structure constant, = 2/3 and 
rut are the charge and mass of the top quark, niz is the Z boson mass, and = 1 — 4 sin^ t^w 
with being the weak-mixing angle. The three terms in Equation ^ correspond to the 
contribution from annihilation via a photon, a Z boson, and their interference. Using the 
published limits of CDF on FCNC, one can derive the following model-dependent limits at 
95% CL I, |]: ^2 < 0.176 and 4 < 0.533. 

^Throughout this paper, charge conjugate states are imphed. 
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In principle, a large FCNC coupling could not only lead to the associated production of 
a top plus a light quark at LEP2, but also to sizable branching ratios of the top quark into 
7c(u) or Zc(u). This analysis uses only the t bW channel. The reduction of the branching 
ratio BR(t — > bW) due to possible FCNC decays is taken into account in the results section. 



2 Data and Monte Carlo Samples 

The present analysis is based on data collected by the OPAL detector [l^ from 1998 to 2000 
at centre-of-mass energies between 189 GeV and 209 GeV. OPAL is a multipurpose high 
energy physics detector incorporating excellent charged and neutral particle detection and 
measurement capabilities. The search presented here uses 600.1 pb~^ of data collected at 
high energies for which the necessary detector components were required to be operational 
while the data were recorded. In addition, 11.3 pb~^ of calibration data were collected at 
a/s ~ mz in 1998-2000 and have been used for fine tuning of the Monte Carlo simulation. In 
this paper, the data sample recorded in 1998 at ^/s ~ 189 GeV is analysed in one bin, while 
the data from 1999 are divided into four samples at y/s ~ 192, 196, 200 and 202 GeV. The 
data collected in 2000 is analysed in two samples of mean centre-of-mass energies of about 
205 and 207 GeV. 

A variety of Monte Carlo samples were generated for the evaluation of the detection 
efficiencies for single top production and SM background processes. In all samples, the 
hadronisation process is simulated with JETSET 7.4 [|ll|] with parameters described in Refer- 
ence |12] and the W boson mass is set to m-yv = 80.33 GeV/c^. For each Monte Carlo sample. 



the detector response to the generated particles is simulated by a GEANT3 based package [|T^ . 

The main Monte Carlo generator used for the description of our signal is PYTHIA [ p| , 
which produces tc(u) via an s-channel exchange of a Z boson. The top quark decays into a 
b quark and a W boson before it can form a bound state or radiate gluons. A colour string 
is formed between the b and c(u) quarks to form a colour singlet. All couplings and quark 
fragmentation parameters for e^e~ tc(u) are set as in Z decays to quark pairs. For an 
evaluation of systematic errors associated with the Monte Carlo modelling, the signal is also 



modelled with a different PYTHIA process and the EXOTIC generator [|1^]. This other PYTHIA 
process is based on a model ||15| for the production of a horizontal gauge boson, called 
R", with the decay R° —>■ tc(u). The EXOTIC generator was developed for pair or single 
production of heavy and excited fermions. Here, the top quark is the heavy fermion and its 
production is associated with a c or u quark. A sequential decay model is assumed with all 
couplings to the known gauge bosons set to the SM expectations. For all three generators, 
samples for three different top quark masses (169, 174, and 179 GeV/c^) have been generated. 
The signal Monte Carlo samples used for the reaction e"'"e~ tc(u) encompass a wide range 
of schemes for the form of the FCNC couplings, the angular distributions of the final state 
particles, and the parton shower parameters. 

The background processes are simulated, with more statistics than the data collected, 
using the following event generators: PYTHIA, KK2F and HERWIG [|17| for (Z/7)*^ qq(7); 



grc4f [|T|], KORALW [|l9l, and EXCALIBUR [|0l for four-fermion (4f) processes; and HERWIG, 



PHOJET [|21|, and Vermaseren for two-photon scattering. 
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3 Event Selection 



The searches for single top events e"'"e~ tc(u) — > bWc(u) are sensitive to the leptonic and 
hadronic decays of the W boson: W iu£ and W qq. The leptonic channel is a clean 
final state with specific topology and kinematics; it is characterised by two hadronic jets, 
one isolated lepton, some missing energy (carried away by the neutrino), and the presence 
of a b-hadron decay. While the hadronic channel is not as clean as the leptonic channel, 
it is statistically significant because BR(W hadron) ^ 68% and BR(W — > iui) ~ 32% 
{£ = e, /i, and r). The hadronic channel is characterised by four hadronic jets with specific 
topology and kinematics, large visible energy, and the presence of a b-hadron decay. Common 
search procedures are applied to both channels. The event selection begins with loose global 
preselection criteria designed to remove most of the two-photon and low multiplicity events. 
To obtain optimal resolution for single top candidates, kinematic fits are performed to reject 
badly reconstructed events and background which are not compatible with the topology 
of single top events. Consequently, the event selection is followed by detailed preselection 
cuts for both the leptonic and the hadronic channels. The final candidate events are then 
identified using relative likelihood functions. Each step will be described briefly in the 
following subsections. 



3.1 Global Event Selection Criteria 

Events are reconstructed from tracks in the central tracking system and energy clusters in 
the electromagnetic and hadron calorimeters, using selection criteria which are the same as 
those used for the OPAL Higgs analysis Because of the presence of jets in a single 



top event, general multi-hadronic preselections are apphed. Each event must qualify as a 
multi-hadronic final state according to the criteria of References p4| , These cuts remove 
events with low multiplicity or little visible energy and reject effectively two-photon and pure 
leptonic events. 

The final state particles and clusters are grouped into jets using the Durham algo- 



rithm |2^. These jets are used as reference jets in the following assignment procedure. 
In calculating the visible energies and momenta, E^s and pVis, of the event and of individual 
jets, corrections are applied to prevent double-counting of the energy of the tracks and their 
associated clusters . 



3.2 Lepton Identification 

Lepton identification for the leptonic channel relies primarily on the isolation criteria of 
a prompt charged particle. Isolated leptons are identified using the Neural Network (NN^) 
described in Reference ||28| . The NN^ uses all tracks in an event with |p| > 2 GeV/ c which are 



considered one-by-one in decreasing order of momentum. They are used as seed tracks and 
all tracks and unassociated clusters within 10° of the seed track define the lepton. Afterward 
the leptons are classified as one-prong or three-prong candidates depending on the number 
of tracks within the 10° cone. Around the seed candidate an annular cone of 30° is drawn 
concentric with and excluding the 10° narrow cone. This serves to define the isolation criteria 
of the lepton candidate. The NN^ provides a distinctive signature for high energy leptons 
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from the particle flow in the annular and narrow cones. However, this procedure is flavour 
blind; the main interest is to retain high identification efficiency. Thus, the NN^ topological 
identification is sensitive to the detection of electrons, muons, and taus with efficiencies of 
84%, 84%, and 75%, respectively. The probability of misidentifying a hadron from a parton 
shower as a lepton is around 1% for NN^ > 0.75. The main source of misidentified leptons 
comes from low-multiplicity gluon jets. 

The lepton with the largest NN^ output in every event is taken to be the lepton of the 
t — s> bW — > b decay. In order to improve the performance of the kinematic fits, a simple 
identification is used to determine the mass (flavour) of the lepton candidate. First, all three- 
prong candidates are classified as taus. Then, a lepton is classified as an electron if Ve > 0.5, 
Ee > 20GeV, and cos9e^u < 0.25, where Ve is the standard OPAL electron identification 
probability Eg is the energy of the lepton, and Oi-u is the opening angle between the 
lepton and the missing momentum vector. Of the remaining leptons, the candidates with 
Ve < 0.5, Ei > 20 GeV, and cosOg^i, < 0.25 are classified as muons, while the others are 
labelled as taus. 



3.3 Event Kinematics 

At a centre-of-mass energy of ^/s ^189 GeV the top quark is produced close to threshold. 
As the top quark is nearly at rest, the W boson in the e"'"e~ — > tc(u) — > bWc(u) reaction has 
almost constant energy i?w — ('^t + '^w ~ '"^b)/2'^t, which also leads to fixed energy for 
the b quark Ei, ~ (m^ — + m^)/2mt. With increasing centre-of-mass energy this unique 
kinematic signature gets diluted. 

These specific kinematic properties are exploited by using kinematic fits. First, the event 
is constrained to pass a 4C kinematic fit which ensures that the energy and momentum are 
conserved]^. The 4C kinematic fit is employed to remove badly reconstructed events and 
events with missing particles along the beam pipe. The probability of the 4C fit is thus 
required to be larger than 10~^. 

To obtain optimal resolution for the reconstructed candidates and performance for the 
jet assignment, we use additional kinematic fits which enforce energy and momentum con- 
servation and impose the appropriate mass constraints. These fits are referred to as the 6C 
fits with the tc(u) or the WW hypothesis: 

• e+e" — > tc(u) — y bWc(u): W boson and top quark invariant mass constraints. 

• e+e" — >-WW: two W boson invariant mass constraints. 



In the 6C kinematic fits the W mass and the top quark mass are fitted with a soft 
constraint, approximating the Breit-Wigner shapes by Gaussian resolution functions. As for 
the 4C fit, we still refer any mass constrained fit as a 6C fit for b£z/^c(u) events. To ensure 
that the kinematic properties of the event candidates match our signal process, we require 
P(6C) > IQ-^ for the tc(u) 6C fit. 

^For semileptonic events, there are three unmeasured variables corresponding to the neutrino momentum 
so that the effective number of constraints in the leptonic mode is one. Nevertheless, any fit which implies 
that energy and momentum are conserved for both the leptonic and hadronic channels will be referred to as 
a 4C fit. 
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3.4 B-Tagging 

The dominant background in this analysis comes from WW events. In e"'"e~ — >WW events, 
the only heavy quark commonly produced is the charm quark. The production of bottom 
quarks is highly suppressed due to the small magnitude of |Vub| and l^cbl and the large mass 
of the top quark. Furthermore, since the top quark is expected to decay into a bW pair, the 
tagging of jets originating from b quarks plays an important role in single top production 
searches. The jet-wise b-tagging algorithm, which has been developed for the Higgs boson 
search, uses three independent b-tagging methods: (1) lifetime tag, (2) high-px lepton tag, 
and (3) jet shape tag. These three methods are combined using an unbinned likelihood 



method to form a single discriminating variable for each jet The b-tag becomes im- 

portant for higher centre-of-mass energies because the kinematic situation changes and the 
signal is less well separated from the WW background. 

3.5 Jet Assignment 

In the hadronic channel, the correct assignment of particles to jets plays an essential role in 
reducing four-jet like backgrounds. There are twelve possible combinations to assign two jets 
to the W boson, the third jet to the b quark, and the fourth jet to the light flavoured quark. 
Therefore a discriminating variable is calculated, which is a combination of the 6C kinematic 
fit probability and the b-tag variable, in order to find the best matching combination to the 
signal hypothesis. The 6C fit helps to identify the two jets coming from the W and to find 
the third jet which matches kinematically to form the invariant top quark mass. In addition 
the b-tag variable helps to identify if this latter jet is a b-jet. The jet assignment which 
yields the largest 'P(b-tag, 6C) is used to choose the jet/quark assignment. 'P(b-tag, 6C) is 
calculated as: 

P(b-tag.6C)- V(.m-t.,) 



V{6C) P(b-tag) + [1 - P(6C)] [1 - P(b-tag)] ' 

where 'P(b-tag) is the b-tag variable and V{6C) is the probability from the tc(u) 6C fit. 

In the leptonic channel, the correct jet/quark assignment plays an important role in 
reducing signal-like background topology. In e+e~ tc(u) — >• b£z/£c(u), there are only two 
possible jet assignments. One of the jets must come from the hadronisation of the b quark 
and the other from the light flavoured quark. The bottom jet is taken to be the one with 
the largest ^(b-tag, 6C) from the tc(u) 6C fit. 

With the jet assignment method described here, a Monte Carlo study shows that the 
rate of correct b-jet (non b-jet) assignment at \fs = 189 GeV is about 96% (94%) and 84% 
(73%) for the leptonic and hadronic channels, respectively. 

3.6 Single Top Candidate Preselection 

To help further reduce the background after the global event selections, the kinematic fits, 
the b-tagging, and the jet assignment, individual preselection criteria are enforced for both 
the leptonic and the hadronic channels. 
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3.6.1 Preselection: Leptonic Channel 

The following preselection cuts are applied in order to reduce background with a different 
topology to our signal process: 



1- ^lopton > 1, where A^iepton is the number of lepton candidates as described in Section ^ 



2. I cos6'miss| < 0.9, where | cos6'miss| is the cosine of the polar angle of the missing mo- 
mentum vector. This cut rejects a large portion of the qq(7) background. 

3. Mvis/i/i > 0.20, where Mvis is the invariant mass calculated from the visible energy 
E^is and the visible momentum p^is of the event. 

4. \Pmiss\/y/s < 0.50, where Pmiss is calculated from the visible momentum {pmiss = —Pvis)- 
This cut reject events with large missing momentum, such as qq(7) background when 
the photon escapes detection. 

5. 0.20 < J2 \pt\/\/s < 0.90, where X) \Pt\ is the scalar sum of the transverse momentum 
components for all the good tracks and unassociated clusters. This cut prevents the 
visible momentum being toward the beam direction and rejects non-radiative qq events 
with no missing energy. 

6. NN^ > 0.75, where NN^ is the primary lepton Neural Network output as described in 
Section I 



After the preselection the background is well described by 4f and qq events. Other final 
states, such as two-photon events, are negligible. The main background (around 95%) is due 
to WW — > qqiue events. The fraction of events with four quarks in the final state selected 
with the leptonic preselection criteria is negligible. 

3.6.2 Preselection: Hadronic Channel 

The following preselection cuts are applied in order to select only four-jet like events: 

1. The event must contain at least 15 charged tracks. 

2. The maximum energy of any electron or muon found in the event (identified as de- 



scribed in Reference ||3^) must be less than 40GeV. 



3. The radiative process e"'"e — > (Z/7)*7^ qq7 is reduced by requiring that the effective 
centre-of-mass energy ^/V' |2^ be at least 150 GeV. 



4. The Durham jet resolution parameter 7/34, at which the number of jets changes from 
three to four, is required to be larger than 0.001. 

5. The (Z/7)*— > qq background is further suppressed by requiring that the event shape 
parameter C [0, which is close to one for spherical events, is larger than 0.4. 

After the preselection the background is well described by 4f and qq events. The expected 
background is composed of 41% (70%) of 4f and 59% (30%) of qq processes at ^/s ~ 189 
(207) GeV. Other final states, such as two-photon events, are negligible. The fraction of 
qqiui events selected by the hadronic preselection criteria is less than 1%. 
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3.7 Likelihood Selection 



The final separation of tlie signal from the background is achieved with a conventional multi- 



variable relative likelihood function ||32|| : 



^ = p Ij" ^ with v = Up.- 

' signal ~r / background j 

The template (or reference) histograms of the input variables, Pi, are used as the probability 
density functions for the calculation of Psignai and Pbackground- We rely on Monte Carlo events 
to compute the probability density functions. 

3.7.1 Likelihood: Leptonic Channel 

For each event satisfying the qqiue preselection cuts, a binned likelihood function is con- 
structed, with one class for the signal and one for the 4f background. The relative likelihood 
is calculated using the following variables: 

Ec{u)'- The energy of the light flavoured jet. 

Mqq'. The invariant mass of the qq system after the 4C fit. 



= y-E^eam ~ iPtu)"^'- Ps^udo mass of the ^Vi system after the 4C fit, calculated from 
the beam energy and the momentum of the pair. 

In^/i2: The logarithm of the Durham jet resolution parameter at which the number of 
reconstructed jets passes from one to two. 

+ M^^: The sum of the di-jet and iu£ invariant masses for the 6C fit under the WW 
hypothesis. 

b-tag: The b-tag variable of the selected bottom jet. 

7-*(b-tag, 6C): The discriminant variable which combines the b-tag variable and the tc(u) 
6C fit probability. 

Jets tagged as light flavoured jets in the background from SM processes have much higher 
values of -Ec{u)- The second and third variables offer good discrimination for WW — *■ qq^i^^ 
events since they exploit the specific angular distribution of signal events. The variable yi2 
exploits the unique jet distribution in tc(u) events. To further remove the background from 
semileptonic WW decays, we use M^^ + Mf^ . Finally, we use the b-tag variable and the 
P(b-tag, 6C) of the b-jet to separate bottom-less events. Figure shows the distributions 
of the input variables for data, the SM background, and the simulated single top signal at 
^/s = 189 GeV. There is good agreement between data and Monte Carlo distributions from 
background processes. 
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3.7.2 Likelihood: Hadronic Channel 

For each event satisfying the qqqq preselection cuts, a binned likelihood function is con- 
structed, with one class for the signal and two for the qqqq and the qq backgrounds. The 
relative likelihood is calculated using the following variables: 

X^iQC fit): The of the tc(u) 6C kinematic fit. 

Ec{u)/ E^is- The ratio of the energy of the c(u) jet and the total visible energy. 
Thrust: The value of the thrust for the event ||11|| . 
b-tag: The b-tag variable of the selected bottom jet. 

cos(Z(pwqi? Pwq2)): The cosine of the angle between the two jets tagged as decay products 
of the W boson. 

The thrust variable exploits the different event topologies between signal and back- 
grounds. The b-tag variable is used as an effective likelihood input because the top quark is 
expected to decay into a b quark. The other three variables exploit the specific kinematics 
of signal events. Figure |^ shows the distributions of the input variables for data, the SM 
background, and the simulated single top signal at ^/s = 189 GeV. There is good agreement 
between data and Monte Carlo distributions from background processes. 

3.7.3 Likelihood: Both Channels 

In Figure ^ the relative likelihood functions for the leptonic and hadronic channels are shown 
for data collected at y/s ~ 205 - 207 GeV and for the SM expectation. No excess of events 
is observed. The likelihood functions for FCNC signal events for an assumed arbitrary 
cross-section of 3pb are also depicted in Figure |[ It can be seen that the expected signal 
contribution is concentrated at high values of C. 

There is no evidence of single top quark production in the data for any y^. Thus, the 
final likelihood cuts are chosen at each value of ^/s so as to minimise the expected upper 
limit on the signal cross-section and thus to maximise the expected exclusion sensitivity. The 
number of selected data and expected SM background events as a function of the centre-of- 
mass energies are shown in Table |I]. 

4 Systematic Errors 

4.1 Signal Efficiencies and SM Backgrounds 

Sources of systematic uncertainties are investigated for their effect on the signal detection 
efficiencies and the SM backgrounds. They are listed in Table |^ for three of the energy bins 
and are discussed below. All checks were performed for all centre-of-mass energies. Possible 
color reconnection and Bose-Einstein effects were not investigated. 

The errors on the background and signal rates from the modelling of the preselection 
variables and of the detector response are a few percent. These uncertainties are evaluated 
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Label 




Lumi. 


Leptonic Channel 


Hadronic Channel 


(GeV) 


(GeV) 


(pb-i) 


Data 


SM Total 


Data 


SM Total 


189 


188.7 


172.1 


3 


4.0 


13 


11.6 


192 


191.6 


28.9 





1.0 


7 


5.1 


196 


195.6 


74.8 


1 


2.9 


6 


6.4 


200 


199.6 


77.2 


3 


2.7 


10 


9.4 


202 


201.6 


36.1 


2 


1.2 


8 


7.5 


205 


205.1 


80.3 


1 


2.0 


11 


10.1 


207 


206.8 


130.8 


6 


3.8 


14 


16.4 



Table 1: The luminosity- weighted mean centre-of-mass energies, the integrated luminosities, 
the number of selected data and expected SM background events at ^/s = 189-209 GeV are 
shown for the leptonic and hadronic channels. 



based on comparisons of the distributions of the variables in the calibration data collected 
at ^/s ~ mz and the Monte Carlo simulation. The effects of detector mis-calibration and 
deficiencies were investigated by varying the jet and lepton energy scales over a reasonable 



range [^. The uncertainties on the energy resolution and the angular resolution were 
also evaluated, but have much smaller effects. A comparison of alternative Monte Carlo 
generators for the background accounts for an additional uncertainty on the background 
rates. The difference between the luminosity-weighted centre-of-mass energies in data and 
the value of ^/s used in the main Monte Carlo samples results in an additional uncertainty 
on the background and signal selection efficiencies due to the use of an energy constraint in 
the kinematic fits. Lepton identification accounts for an extra uncertainty for the leptonic 
channel. 

One of the dominant errors in both analysis channels arises from the b-tagging. Recent 
improvements in the knowledge of heavy quark production processes and decays, such as the 
b-hadron charged decay multiplicity and the gluon splitting rate to heavy quarks, are taken 



into account in the analysis by reweighting Monte Carlo events [^]. The sensitivity to the 
b-vertex reconstruction was assessed by degrading or improving the tracking resolution in 
the Monte Carlo. It was found that changing the track parameter resolutions by ±5% in the 
Monte Carlo simulation covers the range of possible differences between data and simulated 
events. Overall it leads to an uncertainty of 3.8-8.4% for the b-tag rates of background 
and signal events. The finite size of the Monte Carlo samples used in this analysis results 
in an additional uncertainty of a few percent for the background and the signal selection 
efficiencies. 

All the different systematic effects for the background and the signal efficiencies are 
treated as being independent. The total uncertainties on the background and signal rates, 
for both the leptonic and the hadronic channel, are in the same range and show small 
dependencies on the centre-of-mass energy. For each centre-of-mass energy, the systematic 
errors are included in the calculation of the cross-section upper limits. 
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Source of 


Leptonic Channel 


Hadronic Channel 


Systematic Error 


Aefficiency 


Abackground 


Aefficiency 


Abackground 


Preselection 


1.0/1.4/1.2 


2.0/2.2/1.9 


1.0/0.4/0.3 


3.2/1.4/0.5 


Detector Response 


1.0/2.2/3.4 


1.4/1.7/1.2 


0.6/2.0/1.5 


1.0/1.0/3.0 


Background 


-/-/- 


6.8/7.4/6.9 


-/-/- 


5.0/5.0/5.0 


yiin MC 


1.1/2.4/1.9 


2.7/2.1/2.3 


0.6/0.8/0.5 


1.5/1.3/1.4 


Lepton ID 


4.4/5.0/4.8 


3.0/4.0/3.5 


-/-/- 


-/-/- 


b-tagging 


4.2/6.6/5.2 


7.8/5.6/7.0 


3.8/5.3/5.2 


6.9/5.5/8.4 


MC Statistic 


2.2/2.3/2.1 


1.5/2.0/1.6 


2.0/1.8/1.8 


5.4/5.0/4.8 


FCNC modelling 


7.2/8.3/3.5 


-/-/- 


7.9/6.6/5.0 


-/-/- 


Total 


9.8/12.5/9.1 


11.5/10.9/11.0 


9.1/8.9/7.6 


10.7/9.2/11.4 



Table 2: The relative systematic errors (in %) on the signal reconstruction efficiency and 
on the background modelling for ^/s = 189/200/207 GeV. 



4.2 FCNC Modelling 

Several methods of producing FCNC can be compared. A comparison of the results obtained 
with the EXOTIC and the PYTHIA samples described in Section ^ allows an estimate of the 
uncertainty due to the model used for the signal process. The difference is taken as a 
modelling uncertainty on the simulation of signal events. It is summarised in Table ^ 
The main disparities between all the generator schemes are the angular distributions of the 
particles produced in the final state and the parton shower modelling of the initial quarks. 
This latter effect gives rise to one of the largest uncertainties on the signal reconstruction 
efficiency. 

4.3 Top Quark Mass 

The largest systematic uncertainty comes from the sensitivity of the event selection to the 
assumed value of the top quark mass. In the analysis we assume the mass of the top quark 
to be 174GeV/c^. To take this dependency into account, the variation of the reconstruction 
efficiency is investigated using Monte Carlo events with mt= 169 and 179GeV/c^ incor- 
porating the experimental systematic and the FCNC model uncertainties described in the 
previous sections. The dependence of the reconstruction efficiencies on the top quark mass 
for the leptonic and hadronic channels is summarised in Table |^. 



5 Results 

No evidence for single top quark production is observed in e"*"e~ collisions at centre-of-mass 
energies between 189-209 GeV. Limits on the single top cross-section have been derived at 
the 95% CL from the measurements of the number of observed events, the reconstruction 
efficiencies, and the integrated luminosities [^]. The upper limit calculations for each in- 



dividual centre-of-mass energy are summarised in Table 0. Those results include both the 
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(GeV) 



189 
192 
196 
200 
202 
205 
207 



169GeV/c2 



7.5 
7.5 
7.1 
7.1 
6.6 
5.9 
5.8 



10.3 
15.3 
12.8 
14.7 
17.7 
14.4 
12.8 



obs. 
95 



0.30 
0.99 
0.39 
0.55 
1.00 
0.48 
0.47 



nit 



174GeV/c2 



9.1 
9.5 
8.7 
8.0 
7.5 
7.0 
6.7 



12.8 
18.0 
14.7 
16.0 
18.6 
15.7 
15.4 



a, 



obs. 
95 



0.24 

0.81 
0.33 
0.50 
0.93 
0.43 
0.40 



179 GeV/c^ 



6.1 
6.9 
7.2 
7.0 
6.9 
6.2 
6.1 



10.0 
14.9 
12.1 
15.1 
17.3 
13.9 
13.6 



obs. 
95 



0.33 
1.04 
0.40 
0.55 
1.00 
0.49 
0.45 



Table 3: The reconstruction efficiencies for the leptonic (e^) and the hadronic (e^) channels 
are shown. The overall measured 95% CL upper limits on single top production cross- 
section (agg*^') are reported. The statistical and systematic uncertainties are included in the 
calculation of the upper limits. The efficiencies (in %) and the limits on the cross-section (in 
pb) are shown as a function of the centre-of-mass energy for = 169, 174, and 179 GeV/c^. 
These results assume a 100% branching fraction of the top quark into bW. 



statistical and systematic errors and are valid under the assumption that = 169, 174, 
and 179GeV/c2 and that BR(t ^ bW) = 100%. The CDF limits constrain the t ^ Vc{u) 
FCNC branching ratio to be smaller than about 36% |^ for ^ = 7 or Z, so that in a pes- 
simistic scenario the OPAL efficiencies and cross-section limits quoted in Table ^ should be 
rescaled by 64%. 

The combination of all the data can be used to determine limits on the anomalous 
coupling parameters and kz- First, the QCD and the ISR effects which modify the Born- 
level cross-section given in Equation ^ must be considered. The QCD correction is taken 



from Section 3 of Reference while the ISR correction is based on Reference fS^]. Overall, 
the QCD and ISR corrections increase the Born-level cross-section by a constant factor of 
about 1.09 for all centre-of-mass energies and produce only a small distortion to the OPAL 
exclusion region in the — kz plane. 

The limits on the anomalous coupling parameters are obtained with the likelihood ratio 
method described in Reference [^]. Each centre-of-mass energy for the leptonic and the 
hadronic channel has been used as an independent channel. The variation of the selection 
efficiencies for different top masses are taken from Table ^ Taking the statistical and sys- 
tematic errors into account the limit on the anomalous coupling parameters in the — kz 
plane have been derived at the 95% CL. The reduction of the branching ratio BR(t bW) 
due to possible FCNC decays derived at each point in the — kz plane is taken into account 
in this generic FCNC production limit calculation. To compare our results with the limits 
from CDF, exclusion regions for mt = 169, 174, and 179GeV/c^ in the — kz plane were 
obtained. The results are shown in Figure ^. They correspond to upper limits of <0.48 
and Kz <0.41 for a top quark mass of mt = 174GeV/c^, which becomes <0.39 (0.60) and 
Kz <0.34 (0.52) for mt = 169 (179) GeV/c^. These exclusions translate into branching frac- 
tion limits of Br(t Zc) + 5r(t Zu) < 9.7/13.7/20.6 % for mt = 169/174/179 GeV/c^. 
All these results are consistent with recent results from the ALEPH Collaboration 
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6 Summary 



A search for single top quark production via FCNC has been performed with 600.1 pb~^ of 
data collected by OPAL in c+e~ coUision at y/s = 189-209GeV. In total, 85 events were 

selected in the data with a SM expectation of 84.1 events. Limits on single top quark cross- 
sections have been derived at the 95% CL. This leads to model-dependent upper limits 
of <0.48 and kz <0.41 for a top quark mass of nit = 174GeV/c^. The limits become 
<0.39 (0.60) and kz <0.34 (0.52) for rut = 169 (179)GeV/c2. 
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Figure 1: Distributions of the likelihood variables for the leptonic channel at ^/s ^ 189 GeV. 
Comparisons between the data, the SM 4-fermion (light grey), and qq backgrounds (grey) 
are shown. The dashed line represents single top MC events with rrit = 174GeV/c^ and an 
arbitrary cross-section of atop = 3 pb. 
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OPAL 




Figure 2: Distributions of the likelihood variables for the hadronic channel at ^/s ^ 189 GeV. 
Comparisons between the data, the SM 4-fermion (light grey), and qq backgrounds (grey) 
are shown. The dashed line represents single top MC events with mt = 174GeV/c^ and an 
arbitrary cross-section of atop = 3 pb. 
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Figure 3: Distributions of the likelihood variables for (a) the leptonic and (b) the hadronic 
channels. The comparison between the data collected in 2000 at y/s ^ 205 GeV and ^/s ~ 
207 GeV, the SM 4-fermion (light grey), and the qq backgrounds (grey) is shown. The dashed 
line represents single top MC events with rrit = 174GeV/c^ and an arbitrary cross-section 
of (J top = 3pb. 
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OPAL 




Figure 4: The light grey region shows the OPAL exclusion region at 95% CL in the — 
plane for rrit = 174GeV/c^. The exclusion curves for different values of top quark masses 
are also shown. The hatched area shows the CDF exclusion region [^. The OPAL limits 
include QCD and ISR corrections to the Born-level cross-section defined in Equation 
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